.  '  AD-A099  1%%  VIRGINIA  POLYTECHNIC  INST  AND  STATE  UNXV  BLACKSBURG  D— ETC  F/G  12/1 
"  AN  ALGEBRAIC  APPROACH  TO  SUPER-RESOLUTION  ADAPTIVE  ARRAY  PROCES— CTC(U) 

APR  SI  U  A  CADSOtf;  T  p  BRONEZ  N0001A-B0-C-0303 

UNCLASSIFIED  _ _ _  NL 

END 

€-81  I 

OTIC  | 


BIK  RIE  C0PV-  AD  A  0  99  1 4^», 


JECUWTV  CLASSIFICATION  OF  THIS  F*i  '  O. .  •  Znttf+f) 


REPORT  DOCU^N  ,  OH  PAGE 


RIPOFT  NtfMIC 


title  rm  Mini*) 


91  re  READ  INSTRUCTIONS 

AOC _  BEFORE  COMPLETING  FORM 

A.  GOVT  ACCESSION  NOJ  >  CATALOG  NUI.EE* 


»■  TYFE  OF  *EFO*T  t  FEFIOO  COVCFEO 


An  Algebraic  Approach  to  Super-Resolution  I 
Adaptive  Arrav  Processing  ea  / 

- _ _2i  A 

V  7.  AUTMO*<«  .  M  *3| 

C'  James  A.  /Cadsow 
Thomas  P./Bronez  J 

r  FEFFONMING  ORGANIZATION  NAME  *N0  AOOFEM 

Virginia  Tech 

Department  of  Electrical  Engineering 
Blacksburg,  VA  24061 


Manuscript 


S.  FCNFQFmiNG  ONC.  *EPO"T  numEEH 


I.  CONTNACT  ON  SNANT  NUMBERS; 

M 

L  ^  N00014-80-C-0303 vk 


10.  P*oe*AM  ELEMENT.  FNOJEC".  TASK 
AREA  »  »0«K  JNIT  MUMBEFS 


II  CONTFOLLINO  OFFICE  NAME  »NO  AOOFES*  1IET  HEPONT  OAtg 

Office  of  Naval  Research  (Code  436)  '/  firrii 

Statistics  and  Probability  Program  7~ — » -m)>is«»o>  »»nn? 

Arlington,  VA  22217 _  4 

l«.  MONITORING  AGENCY  name  h  AOOHESS/U  alllmtmrtl  !ram  Can  mil  In  t  ClDc  »l  IS.  SECURITY  CLASS.  lot  IAI»  rwani 

^  "  a  /  s  S®*s  :  ** 


i'  17 


»•.  OISTRI0UTIO*  STATE*C*T  r/w •  - 


Unclassified 


lf.vav 


»S«.  OECL  ASSiPtCATiON-  OOWNG«AO»NG 
SCHEDULE 


d*  v 


APPROVED  FOR  PU8LIC  RELEASE:  DISTRIBUTION  UNLIMITED. 


I  »7.  OlSTRl#UTION  STATEMENT  (ot  rl»#  «Aatr«ct  miormd  in  Stock  20,  ft  dJttimt  from  fto+ori) 


»•.  SUPPLEMENTARY  MOTES 


*t  <Ev  «OPOS  (Coruirntn  on  aid*  if  n«c«aaar?  mm  identity  by  block  nimbi) 

Array  Processing,  Algebraic  Characterization,  Wiener  Method, 
Pisarenko  Method,  Maximum  Likelihood  Method 


20.  a**}1  MAC  T  rCdnrfnu*  jn  r«*«fa«  9ldO  it  nkcmaomrr  md  •dmtlty  by  block  numbmt) 

'  In  this  paper,  an  algebraic  characterization  is  made  or  the  problem  of 
resolving  two  or  more  closely  spaced  (in  frequency  vave  number)  plane  waves 
incident  on  a  linear  array.  This  algebraic  characterization  in  turn  suggests 
a  number  of  adaptive  procedures  for  effecting  the  desired  resolution.  One 
of  these  procedures  is  herein  empirically  shown  to  provide  significantly  bettsr 
performance  when  compared  to  ocher  contemporary  procedures  used  in  array  process¬ 
ing  such  as  the  Wiener  filter,  Pisarenko  and  LML  algoriths.  This  includes  both  a 
Spffpr  frpnitpnrv  rp«?n  1  vf n o  panahil  itv  .and  a  faster  _CQHVereeilC6  rate.  w 


oo  1473  edition  o*  i  mov  5S  is  orscle^e 


S  S  -  ‘  -S*  I*  4*  2  50  ' 


Unclassified 

SfCuP»T*  CU  ASH  f  •CATION  Of  *HiS  •Tat 


iyp 


AN  ALGEBRAIC  APPROACH  TO 
SUPER-RESOLUTION  ADAPTIVE  ARRAY  PROCESSING 


James  A.  Cadzov  i  Thomas  ?.  3ronez 


Department  of  Electrical  Engineering 
Virginia  Polytechnic  Institute  and  Scate  University 
Blaclcsburg,  Virginia  21061 


ABSTRACT 

In  this  paper,  an  algebraic  characterization 
is  made  of  the  problem  of  resolving  tvo  or  more 
closely  spaced  (in  frequency  wave  number)  plane 
waves  Incident  on  a  linear  array.  This  algebraic 
characterization  in  turn  suggests  a  number  of 
adaptive  procedures  for  affecting  the  desired  reso¬ 
lution.  One  of  these  procedures  is  herein  empiri¬ 
cally  shown  to  provide  significantly  better 
performance  when  compared  to  other  contemporary 
procedures  used  in  array  processing  such  as  the 
Wiener  filter.  Piaarenko  and  MLM  algorithms. 

This  includes  boch  a  better  frequency  resolving 
capability  and  a  faster  convergence  rata. 


I.  INTRODUCTION 

An  important  array  processing  problem  is  Chat 
of  determining  the  directions  of  propagation  of 
plane  waves  incident  on  a  linear  array  of  uniform¬ 
ly  spaced  sensors  [lj.  Contemporary  speccral 
analysis  has  been  applied  to  this  problem  and  has 
led  to  the  development  of  a  variety  of  processing 
nechoas  thac  are  able  to  resolve  plane  waves  with 
nearly  identical  directions  of  propagation.  These 
oet.ooas  include  the  Wiener  Filter  method  [21.  the 
'iaxlaum  Likelihood  method  [21.  and,  very  recently, 
the  Pisarenko  method  [31-  This  paper  presents  an 
array  processing  approach  based  upon  an  algebraic 
anarocterization  of  the  irrav  processing  problem. 
This  aooroach  is  shown  to  encompass  the  methods 
mentioned  above  as  well  as  auggescing  alternate 
mecnods . 


2’’  d  sin  J. 

a,  -  - -  ,  1  I  It  1  q  ,  C) 

*  \ 

and  the  { n (n) ;  are  uncorrelated  zero  mean  random 
variables  with  variance  j2,  the  'Afe;  are  the  plane 
waves'  complex  amplitudes,  the  are  phase 

angles  dependent  on  the  sampling  instant,  the  [  f'K.- 
are  the  plane  waves'  directions  of  propagation 
relative  to  the  array,  and  V  is  the  common  wave¬ 
length  of  the  plane  waves.  We  assume  thac  the 
are  all  different.  Clearly,  an  estimate  of  the 
spatial  frequencies  directly  yields  an  esti¬ 

mate  of  the  directions  of  propagation  { f ^ ; . 


The  above  set  of  p  instantaneous  measure¬ 
ments  (1)  is  referred  to  a3  a  "snapshot".  To  aid 
the  estimation  of  the  ufc,  we  utilize  a  number  of 
snapshots  taken  sequentially  in  time.  The  irrav 
data  then  has  rhe  fora 


y„(n)  “  n  <n) 
H  o 


r  Ai.jW“'*, 

k-L  * 


0  <_ti^?-l. 
1  «.  m  <  M 


(2) 


where  m  is  the  snapshoc  index  and  M  is  the 
total  number  of  snapshots  used.  In  this  model,  we 
assume  thac  the  phase  angles  ''P|un;  Are  uncorreiaced 
random  variables  uniformly  distributed  on 
This  description  holds  due  to  the  independence  jc 
the  sinusoidal  sources  and  from  the  approximate 
randomness  of  time-sampling  far  below  the  Svquisc 
race. 


It  will  be  convenient  to  reoresenc  toe  given 
data  in  vector  notation.  The  mL-h  snapshot  '3'  will 
be  represented  bv  the  p  •  1  column  vector 


MODEL  OF  THE  ARSAV  DATA 

Let  us  consider  the  model  of  multiple  plane 
waves  incidenc  on  a  linear  array  of  p  sensors 
uniformly  spaced  d  uniC3  apart  in  which  the 
sensor  measurements  are  contaminated  by  additive 
vru.a  noise.  If  there  are  q  olane  waves,  it 
follows  tnac  ac  anv  tarticular  Instant  in  time, 
t'-.a  array  aaca  ycnj,  3  __  n  _  p-1,  has  the  form 


l.v« 


W«  also  define  the  pure  complex  sinusoid  vector  is 


and  the  noise  vector  associated  with  :ne  n-*  snap- 
snot  is 


n-L 


vnere  the  time  wave 


soaclai  trecuencies  are  aiven 


With  the  above  notation,  we  nav  cvr.oactlv  -ep re¬ 
sent  the  snaosnots  :v  tne  oaca  vector 


-  at  ■ _ i _ u _ i _ i 


X- 


equacion 


jyN 

k*l 


L  <  a  <  M. 


(7) 


The  array  daca  ^  is  random  due  co  its 
dependency  on  the  random  phase  angles  (stem-'  and 
che  contaminative  noise  • i  (n)>.  Assuming  chac 
chese  random  variables  are  pairwise  uncorrelated 
and  invariant  with  respect  co  che  snapshoc  index  a. 
ic  follows  chac  each  data  vector  y^  can  be  inter¬ 
preted  as  being  a  windowed  realization  of  a  wide- 
sense  stationary  random  veccor  process.  The  mean 
value  of  this  process  is  the  zero  vector,  while  its 
associated  p  <  ?  covariance  aacrlx  is  specified  by 


s<y  y^.'- 

1  QJ*  O 


3  2l 


,  -,?k  iak 

ic— 1 


(8) 


where  is  che  p  «  p  idenclty  aacrlx  and 

■  &  *  Ati ;  -  is  che  power  of  the  kc!l  plane  wave. 
Since  the  random  veccor  process  is  wlde-sense 
scationary,  the  covariance  aacrlx  R  must  be  posi¬ 
tive  semi-definite,  Toeplicz,  and  Hermltlan.  Ve 
shall  now  give  an  algebraic  approach  co  identifying 
the  plane  wave  frequencies  •’ ■ ,  based  upon  the 
structure  of  the  data  v^  and  che  associated  co- 
variance  matrix  R  . 


li3^  •  In  either  of  chese  cases,  it  is 
intuitively  desirable  co  select  a  coefficient 
veccor  which  is  nearly  orthogonal  to  each  of  the 
data  vectors  in  some  well-defined  manner.  Once 
such  a  coefficient  veccor  has  been  obtained,  che 
plane  wave  frequencies  are  determined  by  examina¬ 
tion  of  the  zeros  of  the  z-cransform  of  this  veccor. 
Specifically,  zeros  that  are  close  co  the  unit 
circle  are  considered  to  be  indications  of  plane 
waves.  Clearly,  closeness  la  a  matter  of  judge¬ 
ment;  it  may  be  conveniently  evaluated  by  search¬ 
ing  for  nulls  in  che  magnitude  of  the  coefficient 
vector's  Fourier  transform  as  given  by 
AU)  -  <a,  sM>  . 

To  obtain  a  mathematical  measure  of  closeness 
co  orthogonality,  ic  is  beneficial  to  introduce  an 
orthogonality  error  vector  *(a)  whose  mch  element 
is  the  inner  product  of  a  with  .  Ue  define 
the  optimum  a  co  be  a  veccor  a*  which  mini¬ 
mizes  some  positive  definite  functional  f  of 
a (a)  .  Hence  we  write 


£(a)  •  [e(l)  ail)  .. 

.  e(M)  ] ' 

where 

«(=>)  *  ‘i.  vffl>  , 

(10) 

and 

f ! e ( a 1 ) ;  *  aln  f  { e  <  a  >  '• 

(11) 

at  A 


III.  .ALGEBRAIC  PROCESSING  APPROACH 

The  approach  co  be  presented  is  dependent  on 
determining  a  nontrivial  p  <  1  vector  a  chat  is 
orthogonal  to  the  noise-free  component  of  each  of 
the  data  vectors  v,,  .  This  orthogonality  is  de¬ 
fined  by  the  general  inner  product  relationship 

1  -  ■  i.  v  -  - 
—  — n  — t, 

*  ‘  A.”a*i5ka  <a.  s  ,  •  ,  1  a  M  .  '9) 

- .  <  —  — k  —  — 

a* . 

imca  me  -14  are  ail  different  and  the  '  s^m-' 
»re  ransom  in  nature,  1  lictie  thought  will  con- 
'.•■.nie  ,'nesei:  that  4  oust  be  orthogonal  to  eacn 
:t  :ne  ;  sinusoid  vectors  9  • 

next  define  t.oe  general  :-transfjrm  Aizl  jt 
m.a  -efficient  vector  i  bv 


-.-ere  d  •  l  z~L  z**  ...  2**m  .  It  is 

then  readilv  snown  tnat  the  irthogonality  of  a 
to  eacn  s_  .  ,  tmolies  tnat  A(zi  must 

lave  ;  finite  zeros  located  in  the  amt  circ.e 

it  •  V.e  p..  ir.ts  CV.  •  4  *'.  1  _  <  _  1  •  '1th  t.nis 
.r.  -inu,  toe  required  smusoida.  frequencies  :an 
le  let er-ir.ed  ?v  examination  >t  me  reros  H  Aici. 


where  A  is  some  prunentiy  chosen  sec  from  which 
che  solution  vector  a*  is  to  be  selected. 

The  inner  product  in  (10)  and  the  functional 
in  (11)  are  general  at  this  point.  We  shall  new 
choosa  in  particular  the  standard  vector  inner 
product  -a ,  y„,'  -  jy, »  and  the, normalized  mean 
square  error  functional  f{e;  *  —  E{  e  ' ;  .  It 
can  be  shown  that  M 

l  ?  0 

f ;  eia)  •  ■  —  E(  ■-a.  y_  ■>  *  *  a  Ji  'll' 

where  R  is  the  covariance  matrix  ’  5) .  The 
functional  (12)  is  co  be  tiin-jiized  according  co 
3s»ae  constraint  such  chac  a*  is  unique  ana  non¬ 
trivial.  Let  us  now  consider  two  possible  con¬ 
straints  . 

(a)  Kvperoiane  Constraint 

The  first  constraint  is  chat  a‘  lies  on  a 
nyperpiane  specified  bv 

A  •  a  ■:  C?:  a  h1**!  i*  •  1.  Ic) 

where  h  is  a  nontrivial  p  «  1  vector.  The 
solution  to  11  with  chis  constraint  can  be  snown 
co  be 

i  •  - . —  R  i  -  -  * 

1  lL  1 

ma  cne  r.miatua  :ri  tenon's  /alee  is  given  b” 


.’on  luea-iistic  ;cnaiti^ns 


■7  r  e  s  e  r.  t 


:*nen..  i:  :  _  t  :r  i :  mere  is  noise 
-r.ere  will  not  -xu:  i  ^efficient  ecccr 
.  ■»  r::vj<onaa  :  ■  .M^n  r  :he  catu  vectors 


f  e 1  * ’ '  *  — 

n  R  *  n 

n )  'uaantic  Constraint 

The  secona  onstramc  is 


ales 


x_. 


I 

I 

i 

i 

i 
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quadratic  surface  specified  by 
A  *  ia  £CP  :  a'  Da  •  U 


(16) 


i 
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M 

/  v  v 

a  -m 

m“l 


(19) 


where  W  is  a  positive  daf inice,  symmetric  p  *  p 
aacrix.  The  solution  to  (11)  with  this  constraint 
can  be  shown  to  be 


j  x '  W  x  , 

.  /  — ain - min  j 


(17) 


It  is  apparent  that  Km  is  unbiased,  Heraitlan, 
but  in  general  not  Toeplit2.  Furthermore,  only 
one  lag  product  from  each  data  vector  is  used  in 
formulating  each  element  of  Ry  .  A  more  desirable 
estimate  is  given  by  the  matrix  jUj  whose  elements 
are 

1^(1. j)  -  c(i-j),  1<1,  j  <_  p  CO) 


and  the  minimum  criterion' 3  value  is 


where  x^ n)  is  the  minimum-eigenvalue  and 

eigenvector  pair  of  . 

These  two  general  solutions  (14)-(18)  encom¬ 
pass  the  three  processing  oechods  noted  in  the 
incroduccion;  (i)  For  Che  choice  h  »  [1  0  ...  0]', 

(1-0  is  the  Wiener  Filcer  solution  [2].  As  in 
linear  prediction,  this  constraint  implies  chat 
the  first  element  of  a’  Is  fixed  and  the  other 
elements  are  unconstrained,  (ii)  For  the  choice 
h  -  s,j,  (lo)  is  the  Maximum  Likelihood  solution  [2J. 
This  constraint  implies  that  A*(z)  has  unity 
gain  at  z  ■  ej“  and  optimally  reduced  gain  else¬ 
where.  (iii)  For  H  "  ip  >  c!le  quadratic  surface 
is  a  hypersphere  of  radius  one,  and  equacion  (17) 
is  a  generalization  of  the  Pisarenko  solution  [31, 

L •* ] .  There  are  several  differences  which  distin¬ 
guish  this  procedure  from  Pisarenko's.  First,  no 
AR2IA  model  is  invoked,  as  is  dona  by  Haykin  [31 . 
Second,  neither  noise  power  removal  nor  matrix 
order  reduction  are  required.  Third,  chis  solution 
is  based  upon  a  minimization  strategy  and  so 
justifies  estimates,  generally  even  non-Toeplitz, 
of  the  covariance  aacrix  R  .  In  the  special  case 
of  u  Toeolitz  estimate,  a  power  Identification 
technique  like  Pisarenko's  can  be  employed,  as  will 
be  shown  later.  Finally,  the  general  constraint 
matrix  W  allows  greater  flexibility  than  does 
the  Pisarenko  method. 


wnere  „  „  , 

1  1  P  i1  1  * 

c(n)  *  w  i  T"r  l  y_d+n)y  U)  ,  0^_n<.p-l 
M  mil  p  n  1»0  a  " 

c(n)  »  c*(-n)  ,  -p  +  l<.n<0  . 

It  is  apparent  that  R^  is  unbiased,  Hermltlan, 
and  Toeplitz.  Furthermore,  it  Incorporates  p-n 
lag  products  in  formulating  the  covariance  element 
c(n).  Therefore  the  variance  ofSjqis  lower  than 
that  of  Sjj  .  Thus,  the  estimate  R»q  is  superior 
to  the  standard  estimate  In  terms  of  its  Toeplitz 
structure  and  lower  variance. 

The  Toeplitz  structure  of  Rj(  has  an  import¬ 
ant  implication  when  used  with  the  hypersphere 
solution.  To  appreciate  this,  consider  a  general 
Toeplitz  Hermician  matrix  with  a  discinct  minimum 
eigenvalue  \mqn  .  An  extension  of  Makhoul's 
findings  [5]  shows  that  the  z-transform  X(z)  of 
the  eigenvector  x  corresponding  to  has  ail 

of  its  zeros  located  on  the  unit  circle.  Ihus  the 
hypersphere  solucion  will  exactly  indicate  the 
presence  of  p-1  plane  waves  if  is  distinct. 

Thus  we  have  a  Pisarenko-like  solucion  and  it  is 
possible  to  apply  a  power  determination  technique 
[4],  [6]  to  separate  the  q  actual  plane  waves 
from  the  p-q-1  spurious  indications  (assuming 
q  <  p). 

Given  an  estimate  of  che  covariance  matrix, 
either  the  hyperplane  or  hypersphere  solutions  can 
be  employed.  We  now  give  simulation  results  for 
these  different  solutions. 


Since  the  Wiener  Filter  solution  has  better 
resolution  than  che  Maximum  Likelihood  solution 
S',  we  snail  hereafter  consider  only  che  hyper- 
oiane  solucion  with  h  »  [i  0  •••  0]'  and  the 
tuauracic  solucion  with  W  •  Ip  (hypersphere 
so lucion i . 

To  summarize  the  development  to  this  point, 
the  algebraic  approach  is  based  on  approximating 
in  orthogonality  condition  between  a  solution 
•■actor  and  each  of  che  daca  vectors.  Thfs 
approach  suggests  many  different  processing 
mechocs.  depending  on  the  choice  of  an  inner  pro- 
aucc.  an  error  functional,  and  a  minimization 
r.'nscraint. 

:v.  zovar lance  MATRIX  estimate 

7 '  -mo  lev  the  hvperolane  and  hypersphere 
solutions  given  lbcve,  an  estimate  of  the  co- 
vir-.jn.e  macrix  .s  required.  A  standard  escimate 


V.  SIMULATION  RESULTS 

To  compare  the  performance  of  chese  processing 
methods,  the  data  vectors  (7)  were  generated  by 
computer  simulation.  The  simulation  modei  corres¬ 
ponded  to  that  chosen  by  Gabriel  [2]  in  his  com¬ 
parative  paper.  Namely,  the  case  of  two  sources 
incident  on  an  array  was  considered .  The  parameter 
selections  were  q  *  2,  p  **  8,  jq-  »  1,  Aq  »  Aa  * 
31.62  (30dB  SNR)  and  3.162  (lOdB  SNR),  rq  •  IS5, 

So  •  22*,  d  *  ’/2,  and  M  »  50  imany  snapsnots)  and 
10  (few  snapshots). 

The  data  vectors  were  analyzed  by  four 
methods:  the  hvperplane  solution  with  estimates 
Rw  and  Rm  ,  and  the  hvpersphere  solution  with  Rw 
ana  Rv  .  3och  the  hvperplane  solution  with  Rv 
and  the  hypersphere  solution  with  R\|  showed  good 
resolution  but  large  spurious  effects.  Results 
for  che  other  two  methods  are  shown  in  Figure  1. 

In  this  Figure,  the  hvperplane  solution  has  been 
evaluated  via  its  Fourier  transform  and  the  hyper- 


sphere  solution  has  been  evaluated  using  the  power 
determination  technique.  Overlayed  solutions  for 
ten  different  realizations  of  the  random  data  are 
shown  to  give  a  sense  of  each  method's  consistency. 

the  results  show  that  both  mechods  work  well 
at  che  high  SNR  with  aany-snapshots.  However,  the 
hyperplane  solution  with  Rjj  performs  very  poorly 
at  low  SNR  with  few  snapshots,  while  che  hyper¬ 
sphere  solution  with  continues  co  give  good 

resolution  and  good  suppression  of  spurious  effects. 
In  general,  che  hypersphere  soluclon  showed  better 
performance  chan  the  hyperplane  solution  over  a 
wide  range  of  conditions. 

VI.  CONCLUSIONS 

We  have  proposed  an  algebraic  processing 
approach  based  upon  approximation  of  a  general 
orthogonality  condition.  This  approach  encompasses 
several  contemporary  high-resolution  analysis 
methods.  One  method  suggested  by  the  algebraic 
approach  has  been  shown  to  provide  significantly 
'oeccer  performance  than  ocher  methods  [2].  Further 


investigation  of  the  algebraic  approach  is  warrant¬ 
ed  in  order  co  fully  exploit  its  potential. 
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